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ABSTRACT 

Most galaxies host central supermassive black holes. As two galaxies merge, the black 
holes also merge. The final single black hole may suffer a kick due to asymmetric 
gravitational radiation and may not be at the centre of the galaxy; off-centre black 
holes may also be produced by other means such as sustained acceleration due to 
asymmetric jet power. We model the main galaxy as a singular isothermal sphere and 
the black hole as an off-centre point lens, and study the critical curves and caustics 
using complex notation. We identify the critical parameters that govern the transitions 
in the topology of critical curves, caustics and pseudo-caustics, and find the number 
of images can be two, three, four and five. We show examples of image configurations, 
including cases where three highly de-magnified images are found close to the centre. 
The perturbation on the image magnification due to the black hole scales linearly with 
its mass in the off-centre case, and quadratically when the black hole is at the centre. 
Such images are difficult to observe unless high-contrast and high-resolution imaging 
facilities (e.g., the Square Kilometer Array in the radio) become available. 

Key words: Gravitational lensing: strong - galaxies: ellipticals and lenticular - galax- 
ies: structure - black hole physics 



1 INTRODUCTION 

Most galaxies host central supermassive black holes (e.g., 
[Giiltekin et al.|[2009| and references therein). As two galax- 
ies merge, the black holes at their centres may also merge. 
The two black holes' orbits first decay through dynamical 
friction when the separation is large. When the separation 
is very small, the black holes can merge efficiently through 
gravitational radiation. However, between these two limits, 
the black holes may stall in their orbital decay. The stalling 
radius is typically at several pc to ^ several tens of pc (see, 
e.g., |Yu||2002| and [Merritt Milosavljevic| [20051 |Colpi fc| 
Dotti|20 Q9'for reviews). 

Whatever brings the two black holes together (e.g., via 
gas processes) , the final single black hole remnant may suffer 
a kick of the order of several thousand km/s due to asym- 
metric gravitational radiation, and thus may be off-centre 
2007). The black hole will oscillate at the 



thors (|Mao et al.||200ll |Chen||2003a|b| [Bowman et al.||2004 



Pretorius 



(e.g.: 

centre of the galaxy, while its amplitude gradually declines 
on the timescale of 1 Gyr ( Gualandris Merritt|20'08 ). An 
off-centre black hole can also be produced by sustained accel- 
eration due to asymmetric jet power ([Tsygan|2007 ) . Such a 
candidate has been reported in M87 (B atcheldor et al.|2010 ) 
with an offset of approximately 12.8 pc. 

Lensing by a single black hole at the centre of a sin- 
gular isothermal ellipsoid has been studied by several au- 



Rusin et al.||2005|. Motivated by observations as discussed 



above, we focus on off-centre black holes and study the criti- 
cal curves and caustics using complex notation. For simplic- 
ity, we model the galaxy as a singular isothermal sphere. 

The outline of the paper is as follows. In ^we rederive 
the basic lensing results for a black hole at the centre. In 
^we study the case with an off-centre black hole, including 
the equations and topologies for critical curves and caustics, 
and then illustrate the image configurations. In ^we discuss 
our results further in connection with observations. 



2 SINGULAR ISOTHERMAL SPHERE PLUS A 
CENTRAL BLACK HOLE 

In this section we investigate a model of a singular isother- 
mal sphere (SIS) plus a black hole at the centre. In complex 
notation ( Witt [1990 ), the lens equation is given by 



/ zz 



l)e" 



(1) 



where z — x ^ iy — re^'^ and z^ — iy^ — r^e^'^'' are the 
complex coordinates of the lens and source plane respec- 
tively and z is the conjugate of z. The term with m is from 
the black hole point lens while the square root term is from 
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the SIS, here m is the black hole mass normalised by the 
total mass enclosed within the Einstein radius. 

For a lensing galaxy with velocity dispersion a = 
250kms~^ at redshift 0.5 and a source at redshift 2, the 
angular Einstein radius is 1.1 arcsec in a cosmology with 
Vim — 0.3, Ha = 0.7 and Hubble constant 70 kms~^Mpc~^. 
The enclosed (cylindrical) mass within the Einstein radius 
is 3.2 X lO^^M©. From the correlation between the black 
hole mass (Mbh) and a ( [Giiltekin et al.||2009J , we find that 
Mbh ~ 2.4 X 10^ M© which implies m ^ l.S x 10~^. The 
scatters in the black hole mass for a given a are quite large, 
0.44 dex in log Mbh- For definiteness we take m = 2.5 x 10~^. 

The above lens equation in polar coordinates can be 
easily transformed into a complex polynomial: 

— rrQe^'^''~'^^ — r — m — (2) 

Due to the circular symmetry, all images must lie on a 
straight line. Without loss of generality we assume all im- 
ages are located on the x-axis and set Lps — ^ (i-e., the source 
is on the positive Xs-axis). Since the solutions for r must be 
real and positive, e*'^ must be restricted to ±1 for this case. 
After some algebra, we find that there are always two im- 
ages, given by 



ri = ^ (^rs + 1 + V'(rs + l)2+4m^ , = 



and 



r2 



^ (l - rs + V(l-rs)2+4m^ 



if : 



(3) 
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The (determinant) of the Jacobian of the lens mapping 
gives the inverse of the magnification: 
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with the derivatives 
dzs . 1 
dz 

and 
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The critical curve is given by J = 0. For our case, in polar 
coordinates, this condition is given by 

{m + r'^){m + r -r'^) =0. (8) 

The first term is always positive while the second term is 
quadratic in r, which yields one real positive solution (the 
negative solution is unphysical) 



(9) 



This critical curve is thus a ring with radius given by the 
above equation. It maps into a degenerate caustic point at 
the origin, as expected from the axis-symmetry. For com- 
pleteness we derive also the magnification of the images (in- 
cluding parity) which is given by 



Ml 



M2 



1 + rs _^ (l + rs)^ + 2m 
2rs 2rsV(l+rs)2+4m' 



(1 



+ 2m 



2^-8 2rsV'(l -rs)2 +4m 



(10) 



(11) 



For m = 0, we recover the familiar magnifications for a SIS: 
/ii = (1 + rs)/rs and fi2 = -(1 - rs)/rs for n ^ 1. In 
astrophysical applications, we expect m <^ 1, so we Taylor 
expand the above expressions into series of m for Ts < 1: 



Ml 



1 + rs 



+ 



(1 + rs 



1 



(l-rs)3 



-.(12) 



The perturbations on the magnifications scale as m^. 

For (fs = and rs > 1, the black hole creates a new 
image which does not exist if we have the SIS alone. The 
approximation for /ii is still valid, but that for the magnifi- 
cation of the new image (/X2) has to be changed to 



M2 



(rs - l)^rs 



(13) 



12 



and 13) for /i2 break down 



Both approximations (eqs. 
when Ts ^ 1. 



3 SIS WITH AN OFF-CENTRE BLACK HOLE 

We now consider the case with a SIS and an off-centre black 
hole. The lens equation is given by 

m z 

Zs — Z — — — — 

z - zo Jzz 



(14) 



where denotes the position of the off-centre black hole 
in complex notation. We choose the coordinate system such 
that the black hole is on the positive x axis with 2:0 = ro > 0. 



3.1 Images 

To solve the lens equation for the image positions it is 
better to switch from the Cartesian coordinates to polar 
coordinates. We again write z — re*'^ and z^ — r^e^"^^ 
(r ^ 0, Ts ^ 0). The lens equation can now be written as 



m 



■ ro 



(15) 



It is interesting to note that this equation differs only by a 
missing r in the last term compared with the binary lens 
equation (assuming one point mass is located at the origin). 

If we clear now the fractions of the equation and take 
the real and imaginary parts of the equation one obtains two 
equations: 



rors sin ifs + rrs s'm{(p - (fs) 
rors cos (fs - rrs cos{(p - (fs) 



= ror sin(^ — ro sin(/?, (16) 
(17) 



= ror cos (p — ro cos cp 
— r + r + m. 



The first equation can be manipulated to yield an expression 
for r: 



rors sin (ps + ro sin (p 
ro sin cp - rs sm{(p - (ps) ' 



(18) 



We may insert this equation into the second equation and 
obtain an equation solely in cp. This equation is equivalent 
to a polynomial of degree 6 in e*'^ and can only be solved 
numerically in general, and so the maximum number of im- 
ages is 6. In practice, we find the number of images can range 
from 2 to 5. For concrete examples, see §3.1| 

For the case when the source position is located on the 
Xg-axis (cps — 0,7r) the equation still factorizes. When the 
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source is on the positive Xs axis (cps = 0), eqs. (16) and (17) 
simplify to 



(f = and r - r(rs + ro + 1) - m + ro + rsro = 0, (19) 



(p — 71 and r +r(rs+ro — 1) 



ro + rsro = 0, (20) 



and 



cos if 



ro 



and 

ro - Ts 

(rprs - m)(ro 
2rors 



+ 



2(ro 



^ 1. 



(21) 



Each equation can have two solutions. However, the latter 
one is bound by the amplitude of coscp and r > 0. In Ap- 
pendix [C] we show in this case, the achievable maximum 
number of images is five. Intuitively, this can be understood 
as follows - there are either zero or two off-axis solutions 
due to symmetry with respect to the x-axis; on the x-axis, 
there are at maximum three solutions: two solutions where 
the black hole and SIS have deflection angles with the same 
sign, and one solution where the black hole and SIS have 
opposite signs. Similar equations can be found for the case 
when the source is on the negative x-axis {cps = vr), but are 
not presented here. 

It is difficult to derive analytically the magnification of 
the images for the general case. However, we can do so when 
the source is located on the (positive) Xs-axis (cps — 0) be- 
cause as we mentioned before the lens equations factorizes 
I 21 ). Appendix [p] gives more details. We men- 



(see eqs. 19 



to 



tion here that the perturbation on the magnification is linear 
with respect to m, in contrast to the m? scaling when the 
black hole is at the centre. 



3.2 Critical Curves and Caustics 

The critical curves are given by the determinant of the Jaco- 
bian J — where J is given by eq. ^ , with the derivatives 

^"^=1-^ (22) 



dz 
and 



Uo 



m \/ zz 



(23) 



In polar coordinates, the determinant of the Jacobian is 
given by 



J= 1 



1 m r 

r r 

where w — 



2rro cos + Tq cos 2(^ 



(24) 



2ror cos + Tq, r > 0, ro > 0. 
The critical curves are given by J = 0, which can always 
be transformed from a two-dimensional problem into a one- 
dimensional one using the parametric representation (see eq. 
8 in 



dz^ 
dz 



Witt|1990| , resulting in the following form 
with ^ a < 27r, 



dz 



(25) 



since the first derivative in eq. ( 22 ) is always real. The equa- 
tion needs to be solved in z for each a in the range to 27r. 
If we switch to polar coordinates we can write the previous 
equation as 



(1 



(ro - re-^'^Y 



2r 



(26) 



If we clear this equation of fractions and take the real and 
imaginary parts of the equation one obtains two equations 
parameterised by a, which give the coordinates in r, cp of the 
critical curve: 



Tq cos(4(/? — a) — 2rro cos(3(/? — a) 
+ (2mr + r^) cos{2(p — a) 
+ (1 — 2r)(ro cos(2(/?) — 2rro cos ip + r^) = 0, 
ro sin(4(^ — a) — 2rro sm{3(p — a) 
+ (2mr + r^) sm{2(p — a) 
+ (1 — 2r)(ro sin(2(^) — 2rro sincp) = 0. 



(27) 



(28) 



For the starting point a = we can obtain analytical results 
for on- axis solutions: 



(p = 0,7T and mr + (1 — r)(ro T r) 



0. 



(29) 



This equation yields one to three solutions on the x-axis. 



However, eqs. (27) and (28) can also have off-axis solutions 
which are more difficult to obtain. Therefore we need to 
disentangle both equations and derive one equation solely in 



r or cos cp. These equations (eqs. Al and A2 ) are presented in 
Appendix [a] Using these starting points one may solve the 
whole critical curve numerically by increasing the parameter 
a from to 27r. 

A transition in the topology of critical curves can take 



place if the following conditions hold (Erdl & Schneider 
19931: 



J(r, (^) = 0, 



a/ 

Or 



0, and 



a/ 

dip 



(30) 



We can use the resultant method ( Erdl &: Schneider||1993 ) 



to eliminate r and e*'^ which yields the condition for the 
transitions of the critical curves. AppendixjB] gives the tech- 
nical details. For m — 2.5 x 10~^, we have 4 transitions at 
ro equals 

(i) ni = 1.26361. 

(ii) rt2 = 0.891231. 

(iii) rt3 = Y^l/4 + m. 

(iv) rtA = Vrn. 

We will illustrate these transitions by examining the shapes 
of critical curves and caustics as we gradually decrease ro in 
a series of figures. As can be seen from Fig.[l] when ro > rti 
there are two disjoint critical curves. One is approximately 
a unit circle associated with the SIS, and the other is a 
small Einstein ring associated with the point lens. These two 
critical curves are mapped into two diamond caustics (see 
the inset). At ro = rti the two critical curves merge into a 
single one (see Fig. |2]), and remains so for m > r > rt2- 
For ro = rt2, the critical curve starts to split into three 
(see Figs. |4] and [5|, with a primary critical curve associated 
with the SIS enclosing two "holes". As the separation fur- 
ther decreases, the two "holes" vanish, leaving behind only 
a single critical curve at ro = rts (Fig. [6|. At even smaller 
separations, the critical curves again split into three sepa- 
rate curves (see Fig.[7|. Another 'peculiar' transition occurs 
when ro = ^/m where the origin becomes part of the caustic 
(see Figs. 8pQ ). This will be discussed in more detail below. 
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3.2.1 Critical curves going through r — 



For a SIS plus a black hole we may have a special transition 
when the critical curve is attached to the origin. In this case 
we have r = and the polar coordinate Lp does not need 
to have a particular value. We investigate eq. (26) for the 
condition for this to occur. 

To do this, we first clear the equation of fractions and 
obtain a polynomial of degree 4 in e*'^ . We may now take the 
complex conjugate of the equation (exchanging Lphy —Lp and 
a by —a) and multiply the new equation by e^*'^. We thus 
obtain two linear independent equations where we can elimi- 
nate e*'^ using the resultant method. Performing these steps, 
we obtain a polynomial in r of the form 4e*^°'r^P22(r) — 0, 
where P22{r) is a polynominal of degree 22 in r (which is 
too cumbersome to present here). The constant term of this 



polynomial is given in the left hand side of eq. (31). 

We note here that r = can not be a generic solution 
of eq. (26). The polynomial p22{r) has a non-trivial solution 
for r = if the constant term in P22 vanishes, i.e.. 



■ 2roe 



+ m = 0. 



(31) 



It is interesting to note that this equation has a valid solution 
for a a rQ ^ m. This means that if the singularity of the 
SIS is inside the Einstein ring of the black hole the critical 
curves are attached to the origin. 

In particular, for the case m — r^ the two inner critical 
curves start to attach to the origin, and is one of the tran- 
sition points described above. In this case we have r = 0, 
a — ^ and Lp — zt7r/2. For ro < \/rn^ a — may no longer be 



, solution for the inner critical curves. However, eq. (31 ) de- 



fines the minimum a value which yields a solution (and the 
starting point) at the origin. That means that the pseudo 
caustic (see the next subsection) becomes part of the caus- 
tics (see Fig. [9]). This is the first case we are aware of in the 
literature where a pseudo caustic merges with a real caustic. 
Furthermore the singularity of the isothermal sphere starts 
to swallow solutions for the parametric representation a. 
Note that for ro ^ the starting point for a moves to 
±7r/2. 



3.2.2 Pseudo- caustic 

When a source crosses a true caustic, the image number 
changes by two. In contrast, when a source crosses a pseudo- 



caustic, the image number changes by one ( Evans & Wilkin- 
son |1998t . The deflection angle for the SIS is not continuous 
at r = due to the singularity, which gives rise to a pseudo- 
caustic. Since z = re*'^, as 



0, the lens equation (14) 



maps into the source plane as 



m 
ro 



^(p <27v. 



(32) 



This is a unit circle with the origin at (m/ro,0), shown as 
the dashed black curve in Fig. [9] 



3.3 Examples of image configurations 

We find that the image number can range from two to five. 



Fig. 10 illustrates the images for 3 source positions (open 
circles, triangles, and diamonds), corresponding 2, 3 and 4 
images. Notice how the pseudo-caustic plays an important 
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Figure 1. critical curves (red) and caustics (black) for a singular 
isothermal sphere (SIS) plus a black hole. The SIS is centred at 
the origin. The black hole is at a distance of ro = 1.35 (indicated 
by a cross). A magnified view of the central region is shown on 
the top right. The left caustic (if further magnified) is similar to 
the one shown in the inset in Fig. |6] 
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Figure 2. Same as Fig. ^except ro = rti = 1.26361. Notice that 
the two isolated critical curves have merged into a single one, so 
have the caustics. 



role. For the case labelled as a square, there are three images 
close to the centre. However, their magnifications are very 
faint, with /x = -3.1 x lO"'^, 4.27 x 10"^ and -8.36 x 10"'^. 
In comparison, the primary image has /x = 1.80. The sep- 
arations between them are of the order of few xlO~^ Ein- 
stein radius (a few milli-arcseconds) for typical galaxy lenses, 
which may be difficult to resolve. 

Fig. [8] shows one five-image configuration for ro — 
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Figure 3. Same as Fig. ^except ro = 1.0. 
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X 

Figure 4. Same as Fig. [l] except ro = rt2 = 0.891231. Notice 
that the critical curve is just creating two "holes" on the right, 
and the caustics are splitting into three segments. 




-1.0 -0.5 0.0 0.5 1.0 1.5 



X 

Figure 5. Same as Fig. [l] except ro = 0.8. The critical curves 
and caustics both have three segments. 
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Figure 6. Same as Fig. ^except ro = rts = 0.5024. There is only 
one critical curve and correspondingly one caustics. 



rt4 = 0.05. In this case, there are 3 faint images created 
by the black hole. The magnifications are —5.98 x 10~^, 
1.16 X 10~^, —9.50 X 10~^ while the brightest primary image 
has /i = 2.06. Notice that the image close to (0.08, —0.05) is 
the negative- parity image associated with the SIS. 



4 DISCUSSION 

In this work, we have studied gravitational lensing by a sin- 
gular isothermal sphere plus an off-centre black hole. We de- 
rived the equations for the images, critical curves and caus- 



tics. We find intriguing critical curves and caustics involving 
the pseudo-caustics. The total number of images for a SIS 
plus a single off-centre black hole can be two, three, four 
or five. In particular, an off-centre black hole can create a 
maximum of three faint images close to the centre (see Fig. 

To the leading order the perturbation on the magnifi- 
cation is quadratic on the primary images if the black hole 
is at the centre and linearly if it is off-centre. Our model is 
simplistic in modelling the primary lens galaxy as a singu- 
lar isothermal sphere. While this appears to be a reasonable 
model for galactic-scale lenses on the scale of few kpcs (e.g., 
Koopmans et al.||2009 ), the central images are likely to be 
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Figure 7. Same as Fig. ^except ro = 0.2. 




Figure 8. Same as Fig. ^except ro = rt4 = 0.05. A source lo- 
cated at (—0.826, 0.449) (green square) has five images (green cir- 
cles): one positive-parity 'primary' image (produced by the SIS) 
at the top left, and four 'central' images. The central images are 
better shown in the inset at the bottom left: the image close to 
(0.08, —0.04) is the negative-parity primary image, while the other 
three are produced by the off-centre black hole. 



very sensitive to the central density profiles (Keeton||2003 



Zhang et al.|2007 ). For a cored isothermal sphere, numerical 
investigations show that the critical curves and caustics re- 
main similar only if the core radius (in units of the Einstein 
radius) is much smaller than m. Magnifications are also af- 
fected. A full investigation of an off-centre black hole in a 
cored isothermal sphere is beyond the scope of this paper. 
An off-centre black hole has been reported in M87 



Figure 9. critical curves (red) and caustics (black) for a singular 
isothermal sphere (SIS) plus a black hole. The SIS is centred at 
the origin. The black hole is at a distance of ro = 0.025 (indicated 
by a cross). A magnified view of the central region is shown on 
the top right. The dashed black curve shows the pseudo-caustic 
corresponding to r = 0, which gives Zs = m/ro — exp(i(/:?) where 
^ < 27r. 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 



Figure 10. Same as Fig. [o] except ro = 0.005. The cross shows 
the point lens position. The open square, diamond and triangle 
symbols indicate three source positions: (xs, ys)= (0.73 1.18), 
(0.39 1.18) and (0.39 0.80). Their corresponding image positions 
are shown as filled symbols. The inset shows the images close to 
to the central region. The red line is the critical curve. There 
are 4, 3, and 2 image positions for the source positions indicated 
by square, triangle and diamond respectfully. The three source 
positions all have an image close to (0.004, 0.0). 
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( Batcheldor et al.|2010| ). The offset is around 12.8 pc, of the 
order of 10"*^ Einstein radius if we put M87 at a typical lens 
redshift (0.5) for a source at redshift 2. The situation will 
be similar to that shown in Fig. [lO] with ro = 0.005. In such 
cases, there is a non- negligible cross-section that multiple 
images at the centre will be formed due to the black hole. 
While these images are rather faint to detect currently, they 
may be observable in the era of the Square Kilometer Ar- 
ray (SKA[^ where its resolution can reach milli-arcseconds 
and the dynamical range can be as high as a million. Since 
the distribution of offsets is unknown, we do not attempt 
a more detailed calculation of the cross-section and proba- 
bilities of seeing multiple central images due to black holes, 
which may also be produced by binary black holes (not yet 
coalescenced) at the centre of galaxies ( |Li et al.||2011 ). 



ACKNOWLEDGMENTS 

We thank Jin An and an anonymous referee for very help- 
ful comments and criticisms that improved the paper. We 
acknowledge the Chinese Academy of Sciences for financial 
support and the hospitalities of the Aspen Center for Physics 
where this work was completed. 



1 www.ska.org 



8 Mao & Witt 



APPENDIX A: THE STARTING POINT OF THE PARAMETRIC REPRESENTATION 



Using eq. (26) and its complex conjugate for a = one can derive one equation in r and another in cosLp using the resultant 
method ( Erdl Schneider]|1993 ). For r we obtain the following polynomial: 

4r^ — 4r'^ + (1 + 4m — 4ro)r^ — 4mr^ + (m + ro)r^ + 4ro(m + ro)r^ — ro(2m + ro)r^ + ro(m — Tq) = 0. (Al) 

For u = cos we obtain 

64roU^ - 32rW - 4ro(32ro - 4m - - 16ro(2ro - 2m - 3)u^ + 4ro(22ro - 11m - l)u^ + 
4ro(m + 6m^ — 4ro — lOmri + 12ro)i^^ + (4ro — 12mro — 24ro + 12m^ro + 16mro — 4m^ — m^)u' 



-4ro{m — 4ro)(m — rQ)u — 4ro(m — Tq)^ 



0. 



(A2) 



Again one sees that when ro = ^/m? r = becomes a solution of eq. ( Al ), which signals one of the transitions in the topology 
of critical curves (see ^3.2.1). 



APPENDIX B: TOPOLOGICAL CHANGES IN THE CRITICAL CURVES 



Following Erdl & Schneider ( 1993| ), the topology of critical curves changes when the following conditions are satisfied 



J = 0, 



0(p or 



0, u = cos(p, w = — 2ror cos (p -\- Vq . 



where J is given by eq. (24). For the derivative with respect to (p, we have 

dJ _ 2mro sin cp (2mr^ + — 3rro + 2rQU 
dip rw^ 

Similarly we find 

_ -m^r + (-1 + r)w^ - m{w - 2rg(l - u^)) _ 
J — - — u 

and 

dJ_ _ 4m'^r^(r - rpu) + + m(3r^ - lOr^rpu + 6rr^u{l - 2v?) + ro(-l + 2u^) + 6rVg(-l + ?>v?)) 
dr r'^w^ 



Eq. ( B2 ) is satisfied for either 



sin Lp — ^ 
or 



u = -(2mr^ - 3rro)/(2rg). 



(Bl) 
(B2) 

(B3) 

(B4) 

(B5) 
(B6) 



For the condition in eq. (B5), we have either (p — ^ or ip — -k. For each case, the conditions J = and dJ/dr — give 
two equations in terms of r, we can use the resultant method ( Erdl Schneider|1993 ) to eliminate r to find the condition for 
topological changes in the critical curve. For ip — u = cos ip = 1, we have 

+ 4mVg(m + ro)(4m + (1 - 2ro)^)^ [m + 4m^ - 4ro - 20mro + 12ro - Smrl - 12rg + Atq] = 0, (B7) 

All the terms in front of the bracket [] are positive - only the bracket term may yield physical solutions, which can be solved 
analytically since it is a quartic equation in terms of ro. Similarly, for = tt we have 



■ 4mVo(m + ro)(4m + (1 + 2roff [m + 4m^ + 4ro + 20mro + 12ro - 8mro + 12ro + 4ro] = 0. 



(B8) 



For physical situations, the black hole mass most likely satisfies m < 1, and all terms in this equation are positive (using 
12ro > 8mro), and so there are no physical solutions. 

For the condition in eq. (B6), the resultant method gives 



32m (1 + 4m - Arlf (m - rg)^(8m^ + 91m^ + 344m^ + 432m^ + (IGm^ + 318m2 + 1659m^ + 2976m^ 
+1296m^)rg + (256m + 2424m2 + 7333m2 + 7248m^ + 864m^)r^ + (1024 + 6208m + 8517m2 - 2272m2 - 
864m^)rg - (3072 + 10304m + 744m2 + 1296m^)r§ + (3072 - 3072m - 432m^)r^° - 1024r^2) = 0. 

This equation has at least two analy tical so lutions from the first two terms in brackets (given below as rta and vtA) 
For m = 2.5 x 10~^, using eqs. ( B7|B9 ) we find four positive physical solutions of ro 



(B9) 



(i) Tti — 1.26361 from the condition in eq. (B7) 



(ii) rt2 — 0.891231, another solution from eq. (B9) in adddition r^s and VtA below 

(iii) rt3 v'l/4 + m. 

(iv) rt4 
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A fifth positive solution ro = 0.00062461 (from eq. |B7| ) does not any give positive solution of r, and so is discarded. These 
four transitions are illustrated in Figures l2] [4] [6] and [8] 



APPENDIX C: MAXIMUM NUMBER OF IMAGES 



To get a maximum of six images, each equation from (19) to (21) must yield two solutions. Since r ^ and —1 ^ coscp ^ 1, 
to have two solutions, eq. (21) must satisfy 



ro > rs ^ 0, m ^ rsro 



(1 + rs 



■ro) 



ro 



Furthermore for eq. (20) to have two positive solutions, we must have the coefficient for the linear term to be negative 
rs + ro < 1 

and the constant term to be positive 
- m - ro + Tsro > 0. 



However, combining the conditions in eqs. (CI) and (C2) we find that the left hand side of eq. (C3) satisfies 



■ m - ro + rsro ^ -rsro 



(1 + rs 



ro 



ro 



■ ro + rsro = (-3 + ro + rs)r 



ro 



ro 



< 0, 



(CI) 

(C2) 
(C3) 

(C4) 



in direct contradiction with the requirement in eq. (C3). In other words, this equation cannot have two solutions, and so the 



maximum number of images is at most five, as we argued intuitively in §3.3| 



APPENDIX D: MAGNIFICATION OF THE IMAGES WHEN THE SOURCE IS ON-AXIS 



When the source is located on the (positive) Xs-axis (cps — ), the lens equation still factorizes (see eqs. 19 to 21); the image 
positions can be derived analytically and we can use eq. (24) to obtain the magnification ji — . 

For = we can find the solution for r by solving the quadratic equation in eq. ( 19 ) and then obtain the magnification 
using eq. (B3) {u = coscp = 1) 



r(r - ro)^ 



_ j-i _ 

—m?r — m{r — roY + (~1 + r)(r — ro)^ 

Similarly for = tt, we find 

r(r + ro)^ 



-m?r — m{r + ro)^ + ( — 1 + r)(r + ro)^ ' 



(Dl) 



(D2) 



where the solutions for r can be found from eq. (20). For m ^ 1, one can Taylor expand these expressions, and find that the 



magnification has a linear perturbation term with respect to the black hole mass (m) for the two outer (primary) images: 



rs±l , 



ro 



Kl-ro+rs)2 



m + 0{m ), 



(D3) 



where the + and — signs are for Lp — and tt (eqs. Dl and D2) respectively. There is one new on- axis image created by the 
black hole whose magnification scales as m^, but we do not give the expansion here. 

The magnifications for the two off-axis images are identical due to symmetry and are given by 

2mro(ro - r^) 



where 



A1A2 



Ai,2 = m(ro - rs)^ - rors(ro -rsT 1)^- 



(D4) 



(D5) 



It is interesting to note that the magnification is of order m. These expressions are only valid when off-axis images exist, i.e., 
when the source is inside a tiny caustic close to the centre. 
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